X-ray natural dichroism and chiral order in underdoped cuprates 
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The origin of the Kerr rotation observed in the pseudogap phase of cuprates has been the subject 
of much speculation. Recently, it has been proposed that this rotation might be due to chiral charge 
order. Here, I investigate whether such order can be observed by x-ray natural circular dichroism 
(XNCD). Several types of charge order were considered, and they can give rise to an XNCD signal 
depending on the stacking of the order along the c-axis. 
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Since its observation over twenty years ago, the na- 
ture of the pseudogap phase i n un derdoped cuprates has 
been a subject of much debate^. Time reversal symme- 
try breaking was proposed by Varma^, and based on this 
prediction, angle resolved photoemission (ARPES) ex- 
periments were done by Kaminski et al™ on underdoped 
Bi 2 Sr 2 CaCu 2 08+5 (Bi2212) that did indeed observe a 
dichroism signal by using circularly polarized light. Their 
signal was consistent with an order parameter like tem- 
perature evolution that set in below the pseudogap tem- 
perature T*. These results, though, have yet to be con- 
firmed by other groups^. Moreover, even if confirmed, 
such a signal could instead be due to chiral order or cer- 
tain types of structural distortions rather than time re- 
versal breaking. On the other hand, subsequent elastic 
neutron scattering experiments did observe time reversal 
breaking in underdoped YBa2Cu30g+a; (YBCO)SJ that 
sets in below T* . This has now been seen in a variety of 
cuprates 7 . Interestingly, this symmetry breaking has not 
been seen by //SRP. 

After the neutron results, a novel Sagnac interferom- 
etry experiment was done by Kapitulnik's groupP that 
observed a Kerr rotation that sets in at a tempera- 
ture below that where the elastic magnetic signal is first 
seen by neutrons. In Bi2Sr2Cu06+<5 (Bi2201), a simi- 
lar signal was seen that set in at the same tempera- 
ture ARPES measurements saw an energy gap develop^. 
Originally, it was suggested that this signal was a signa- 
ture of a tiny ferromagnetic moment, but as Orenstein 
has pointed outP^, a Kerr rotation was seen in the anti- 
ferromagnet Cr2 03 where a structural distortion accom- 
panies the magnetic ordep^l. This occurs since although 

< P > and < M > vanish because of the staggered order 
(where P is the polarization and M the magnetization), 

< (P ■ M)j_ > is non-zero where _L denotes the compo- 
nents of the scalar product perpendicular to the propaga- 
tion vector of the light. That is, the structural distortion 
flips sign with the moment, and so the scalar product is 
invariant to the staggering. A different magneto-chiral 
phase has recently been advocated by Varma and col- 
laborators based on orbital currents to explain the Kerr 
effect in the cuprates 13 . But the idea of chiral charge or- 
der has also been propose d 14 * 15 ^. One motivation for the 
latter proposal is that the onset of the Kerr rotation in 



YBCO does not occur at T*, but rather at a lower tem- 
perature where charg e order has been seen to develop by 
x-ray scatterin g ^ 16 * 17 !. 

To further explore this, it would be useful to identify 
other probes that could detect such chiral order. XNCD 
is a natural one to propose^. This is the difference in 
absorption of left and right circularly polarized light from 
time even processes. It is sensitive to inversion symme- 
try breaking, and as it is site specific, it can be used to 
gain information on the spatial nature of the order. In- 
terestingly, an XNCD signal was seen at the Cu K edge 
in Bi2213i^ that has an order parameter like evolution 
with temperature that matched the ARPES dichroism 
signal. Though many space group refinements of Bi2212 
break inversion symmetry, because of the presence of 
glide planes, the XNCD signal due to the crystal struc- 
ture is zero when the light is directed along the c-axis^Sl 
as in the experiment. A lower symmetry than simply 
the crystal structure would be required to generate an 
XNCD signal, and chiral order could indeed satisfy this 
condition. 

In the present paper, I investigate whether such chiral 
order could be the source of an XNCD signal, and find 
that if the order is such as to violate the two glide planes 
present in Bi2212 (one perpendicular to the Cu02 planes, 
the other parallel), then an XNCD signal does occur. 

I start with a discussion of the Bi2212 space group. 
The basic unit cell is face centered orthorhombic. For 
purposes here, the superstructure modulation observed 
in Bi2212 is ignored, since this involves a translation 
operation that is not relevant to the symmetry consid- 
erations allowing for XNCD. In our previous study^Sl, 
several space group refinements were considered, but for 
the purposes of this paper, two will be studied. The 
first one (Bb2b) of Gladyshevskii and FHBdger^ is non- 
centrosymmetric, the second one (Bbmb) of Miles et alW* 
is centrosymmetric. All refinements have in common the 
existence of two glide planes (with the glides along the 
supermodulation b-axis), one perpendicular to the a-axis 
which runs through the planar oxygen iona^, the other 
perpendicular to the c-axis (midway between successive 
Cu0 2 bilayer units, which is also midway between the 
two BiO planes) . One can easily see that any given order 
must violate both glide planes in order to give a non-zero 
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FIG. 1: Oxygen sites in the four CuC>2 planes of the unit cell 
of Bi2212 (shown are four oxygens per plane). The numbers 
refer to the space group operations in the text, with the two 
different planar oxygen sites shown as unprimed and primed, 
x denotes the location of the coppers. The a-axis runs along 
1-1', the b-axis along 1-3, with a glide plane along 1-3 as well. 
The other glide plane (_Lc) is shown as the dashed line. 



signal, since XNCD involves a sum over all the atoms in 
the crystal associated with the x-ray edge studied. 

XNCD originates from interference between dipole 
(El) and quadrupole (E2) contributions to x-ray 
absorption 18 . This requires inversion symmetry to be 
broken. This contribution is a time even pseudo- 
deviatop22 which for most (but not all) symmetries 
would give rise to optical dichroism as well (a pseudo- 
scalarpSl. There can be a contribution to XNCD as 
well from interference between electric (El) and magnetic 
(Ml) dipoles, but this contribution is very small, which 
was verified in the present case by numerical simulations. 

Numerical results are generated using the FDMNES 
prograrrPSEZl The simulations were done using lo- 
cal density atomic potentials (Hedin-Lunqvist exchange- 
correlation function) in a muffin tin approximation that 
considers multipl e scatt ering of the photoelectron around 
the absorbing sit e! 28 * 29 !. The cluster radius is limited by 
the photoelectron escape deptlP2. Here, fixed cluster 
radii of 3.lA and 4.9A are considered as in Ref. |5D1 but 
in one case, larger clusters of radii 5.5A and 6.5A were 
studied. For all calculations shown, the incident light is 
along the c-axis. 

To proceed, discussing symmetries of the Cu02 planes 
for Bi2212 is in order (Fig. 1). There are two different 
planar oxygen sites, each of which forms rows along the 
b-axis, with the different rows alternating as one moves 
along the a-axis. Associated with this is a glide plane or- 
thogonal to the a-axis. There are also four CuC>2 planes 
per unit cell, composing an upper bilayer and a lower 
bilayer, with another glide plane in between the two bi- 



FIG. 2: (Color online) Three charge patterns considered in the 
text - (a) pattern 1, (b) pattern 2, and (c) pattern 3. Circles 
and squares are oxygen ions, crosses are coppers. Positions of 
the ions are as in Fig. 1. In the calculations, squares have a 
charge excess of 0.1, circles a charge deficiency of 0.1. 



layers. The two oxygen rows are interchanged as one 
moves from each layer to the next, with the upper plane 
of the upper bilayer and the lower plane of the lower bi- 
layer having the same pattern, and the other two planes 
staggered relative to this by half a lattice constant along 
the a-axis. Note that as a consequence of these two dif- 
ferent planar oxygen sites, even in centrosymmetric re- 
finements of Bi2212, inversion symmetry is locally broken 
about the copper sites. As for the coppers, they stagger 
differently, translating along the a-axis by half a lattice 
constant as one goes from one bilayer to the next. This 
can be understood from the eight space group operations, 
which are (1) the identity, (2) two-fold rotation (-x,y,-z), 
(3) first glide (-x,y+l/2,z), (4) second glide (x,y+l/2,-z), 
and then (5-8) the above operations followed by the face 
centering translation (1/2,0,1/2). These operations are 
used to label the oxygen ions in Fig. 1. 

I will now consider charge ordering patterns on the 
oxygens that do not break translational symmetry, since 
this is what the neutron data indicate^, and as mentioned 
above, the translational part of any order is not relevant 
for the symmetry considerations allowing for an XNCD 
signal. Fig. 2a shows a pattern that is the charge ana- 
logue of the oxygen spin pattern that was considered in 
Ref.|H]as an alternate to orbital currents, referred to here 
as pattern 1. Note that this pattern breaks the first 
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FIG. 3: (Color online) XNCD for the various charge patterns 
described in the text at the Cu K edge in Bi2212 (with k 
along the c-axis) for the crystal refinement of Ref. 1211 Also 
shown is the absorption signal itself. The left column is for a 
cluster radius of 3.lA, the right one for 4.9A. The top row is 
for pattern 1, the middle row for pattern 2, and the bottom 
row for pattern 3. Note the differing scales on the right y-axis 
for each row. 



glide. Another pattern, referred to as pattern 2, flips 
the charges in the second row of oxygen ions (Fig. 2b). 
This pattern also violates the first glide. To violate the 
second glide, these two patterns must have the correct 
staggering along the c-axis. For pattern 1, the staggering 
is (-,+,+,-) where the order refers to the plane sequence 
(upper plane, upper bilayer; lower plane, upper bilayer; 
upper plane, lower bilayer; lower plane, lower bilayer) as 
in Fig. 1, with the sign referring to that of the oxygen 
site in the lower left corner of each plane (6', 1,2, 5'). All 
other stackings with an even number of + and - signs 
give a vanishing signal. For pattern 2, a non-zero signal 
under the same condition of an even number of + and 
- signs requires that all planes are in phase (+,+,+,+). 
In both cases, the space group reduces to B121 with four 
operations: the identity, two-fold rotation (-x,y,-z), and 
the face centering translation (1/2,0,1/2) of these two. 

Next, a true chiral pattern is considered based on the 
suggestion of Ref. HU This pattern is the same as in 
Fig. 2b, except it is rotated around the center of the 
four oxygen sites in each plane as one goes from plane 
to plane, with a rotation sequence of (-90° ,0°, 90°, 180°). 



FIG. 4: (Color online) XNCD for the various charge patterns 
described in the text at the Cu K edge in Bi2212 (with k 
along the c-axis) for the crystal refinement of Ref. 1221 Same 
notation as in Fig. 3. 



This pattern, though, preserves the symmetry operation 
(-x+l/2,y+l/2,z+l/2) and thus has a vanishing XNCD 
signal. An alternate pattern presented in Ref. \T5\ has 
the rotation sequence (-90 o ,0°,-90 o ,0°). This pattern 
(Fig. 2c), denoted as pattern 3, breaks all symmetries 
(PI space group), allowing for XNCD. 

Fig. 3 shows the calculated XNCD signal at the Cu K 
edge for the three patterns with a cluster radius of 3.1 A (a 
copper and its five surrounding oxygen ions) in the left 
column, and a cluster radius of 4.9A (which contains 37 
atoms) in the right column. Spin-orbit was included, but 
this only had a modest effect. A charge imbalance of ±0.1 
e~ was assumed on the oxygen sites. Note that flipping 
the charge pattern flips the dichroism signal. Complex 
energy profiles are found. They differ between the various 
charge patterns, and are sensitive to the cluster radius. 
Note that the dichroism in pattern 3 is half that of pat- 
tern 2. This can be understood from the symmetry of 
the two patterns shown in Fig. 2. The magnitude of the 
dichroism is significant, with values up to 4 x 10~ 3 of the 
absorption maximum, which is not only larger than that 
claimed in Ref. Qj|l but well within the range of detection 
of modern x-ray sources. The largest signal typically is 
seen at the pre-edge peak at 8.972 keV, expected since 
this corresponds to Is — 3d excitations (as opposed to 
the edge itself, which corresponds to Is — Ap excitations). 



This is in contrast to the experimental result, which is 
simply a positive peak at 8.99 keV followed by a nega- 
tive peak at 8.998 keV. As discussed in our earlier worlPSl, 
the observed signal seems most consistent with a small 
energy shift between left and right circularly polarized 
light, rather than intrinsic dichroism. Regardless, the 
XNCD energy profile calculated here is very sensitive to 
the nature of the charge order, thus providing a unique 
fingerprint for this order. 

In Fig. 4, results are shown for the centrosymmetric 
refinement of Miles et aZp2l. Generally, the dichroism is 
about an order of magnitude smaller, not surprising since 
this refinement preserves inversion symmetry. Larger 
cluster radii were also run for pattern 2 (5.5A and 6.5A). 
The 5.5A result was similar to the 4.9A one, whereas 
the 6.5A one had a pre-edge XNCD signal which was 



partially suppressed. Still, the XNCD signal is generally 
large enough that it should be detectable if present. 

In summary, I find that XNCD should be an exquisite 
probe of chiral order in cuprates, assuming that the order 
is such to allow a non-zero XNCD signal. Each charge 
pattern calculated has a unique energy profile that should 
be detectable at modern x-ray sources. In fact, a vari- 
ety of potential x-ray measurements could be done that 
would allow for a thorough interrogation of such order, 
including its potential magneto-chiral nature^^U, and it 
would be quite interesting if such experiments were at- 
tempted on a variety of cuprates. 
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